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Abstract 

We investigate a model that the Yukawa coupling form is constructed by two 

kinds of matrix (Mo and Mi). For example, in the 50(10) GUT model, Mq 

and Mi are Yukawa couplings generated by the 10 and 126 Higgs scalars. We 

study how this model can give the observed mass and mixings of quarks and 

leptons. Parameter fitting is fully scanned by assuming all the input data to 

be normally distributed around the center value. 
PACS number(s): 12.15.Ff, 12.10.Kt, 14.60.Pq 
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I. INTRODUCTION 



The grand unification theory (GUT) is very attractive as an unified description of the 
fundamental forces in the nature. However, in order to reproduce the observed quark and 
charged-lepton masses and mixings, a lot of Yukawa couplings are usually brought into the 
model. We think that the nature is simple. So it is the very crucial problem to know the 
minimum number of Yukawa couplings which can give the observed fermion mass spectra 
and mixings. However, if the quark and lepton Yukawa couplings are composed by only one 
matrix 

M u = c u M , M d = c d M , M e = c e M , (1.1) 

the CKM Matrix must be diagonalized, these model is obviously ruled out for the description 
of realistic quark and lepton mass spectra. Therefore, at the unification scale /x = Ax, we 
assume the Yukawa coupling of up quark, down quark and charged-lepton (M°, M° and 
M° ) are composed by two matrices, 

M° = c /0 M ° + c fl Ml (f = u, d, e) (1.2) 

Here, c/ and c/i are real numbers which can be associated with the vacuum expectation 
values (VEV). For example, in the SO(10) GUT model with one 10 and one 126 Higgs 
scalars, the Yukawa couplings of quarks and charged leptons are expressed in the following 
forms [1] [2]: 

M° = c M + dMf, M° = M ° + M°, M e ° = M ° - 3M°. (1.3) 

where M° and M° are symmetric Yukawa couplings. In the previous paper [1], eliminating 
M° and Mi from Eq.(1.2), we obtain the relation 

Ml = c u M° u + c d M° d = c u (M° u + KM° d ), (1.4) 

where 
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CuOCel — C e ()C u i C e QC d \ — C d ()C e i , _•, 

c d = and c u = . (1.5) 

C u oQl ~~ Qo c ul c uO c dl — c d0Cul 

These relations are realized at the GUT scale, but each value of the Yukawa couplings is 
given by the experiment at the weak scale \i = mz- Therefore, we must investigate how the 
mass ratios and CKM matrix parameters change from \x = Ax down to fi = mz- [3] In this 
paper, we distinguish between the values at fi — Ax and fi = mz by using the superscript 
"0" or not. 



II. NUMERICAL STUDY 

Because M°, M d , and M e ° are symmetric at the unification scale fi = A x in the model 
with one 10 and one 126 Higgs scalars, they are diagonalized by unitary matrices U®, U d , 
and C/g, respectively, as 

U°JM° U U° U * = D° u , V*M°V? = D° d and U°J M° e U? = D° e , (2.1) 

where D®, D d , and D® are diagonal matrices which are given by 



, l + tan/3 2 o o\ n o l + tan/5 2 , „ (l 

K = j diag(m^, m u c , m" t ) , D d = - diag(m^, m u s , m v b ), 

v v 



^ e °- 1± ^diagK,mJ,m°), (2.2) 

v 

Here, v(— 174 GeV) is VEV of Higgs , and it is divided into up and down quark (neutrino 
and charged lepton) in the ratio tan/5. Using the Cabibbo-Kobayashi-Maskawa (CKM) 
matrix V q which is expressed as V q = U^U d , the relation (1.4) is rewritten as follows: 

(U^JD^KY = c u Dl + c d V q °D d V q ° T = c u {Dl + <A^ or ). (2.3) 

We take a basis on which the up-quark Yukawa coupling is diagonal in order to compare 
with the experiment values and obtain the independent two equations: 

(m°m°) 2 + (™>°) 2 + Km°) 2 ) 2 [Tr{ ^ {k)}] 2 



A(k) 



m° e ) 



f + K) 2 + k) 2 ) 2 ^ ^» 2 - Tr fa ^ 2 ) 



(m>X) 2 [Tr{^(/,)}] 3 



m°) 



Here H q (n) is the following hermite matrix which is defined by the Yukawa couplings of 
quark: 

H 9 (k) = (D° u + KV°D° d V°W + KV°D° d V ^. (2.5) 

If we find the k which sets A(k) and B(k) to 1 simultaneously, the three Yukawa couplings 
M°, M° and M° can be unified into two matrices. However we don't know precisely how 
to determine these data, especially quark masses. And above procedures depend on these 
ambiguities. So in this paper, we substitute the random numbers which becomes following 
normal distributions [4]: 

\m u (2GeV) | = 2.9 ± 0.6MeV, \m d (2GeV) | = 5.2 ± 0.9MeV, (2.6) 
|m s (2GeV)| = 80 - 155MeV, \m c (m c )\ = 1.0 - 1.4GeV, (2.7) 
m b (m b ) = 4.0 - 4.5GeV, mf rcct = 174.3 ± 5.1GeV, (2.8) 



m pole 



m pole 



= 0.510998902 ± 0.000000021MeV, (2.9) 
= 105.658357 ± 0.00005, m p T olc = 1776.99 ± 0.29MeV (2.10) 



sin 12 = 0.2229 ±0.0022, sin 9 23 = 0.0412 ± 0.0020, (2.11) 
sin 13 = 0.0036 ±0.0007, 5 = (59 ± 13)° (2.12) 

for each mass and CKM mixing parameter 10,000 times. And we estimate the evolution 
effect about the values in Eqs. (2.6) - (2.12) from fi = mz to fi = Ax by using of RGE. [3] 
In this work, we suppose MSSM for tan/3 = 10. Without loss of generality, we can make 
the masses of third generation positive real number. Although the remaining masses are 
complex under the ordinary circumstances, we assume that all masses are real in order to 
simplify the problem. Therefore, there are 16 combinations of the signs of the masses as 
shown in table I. 

As shown in Fig.l, we scan the range A(k) = 1 by changing Im(ft) from -100 to 100 at 2000 
equal intervals. Moreover, we get the maximum and minimum of B(k) on the line of A(k) = 1 



by changing Im(/t) at 5000 equal intervals. Because B(k) is continuous, there is the n which 
sets A(k) and B(k) to 1 simultaneously when Min(£>(/t)) < 1 < Max(S(/c)) as explained in 
Fig. 2. In this way, we draw the histograms in Figs 3,4,5 and 6 which show the distribution 
of input values conforming to the requirements Mm(B(n)) < 1 < Max(S(/c)). The each 
summation of the conforming case is tabulated in Table 2 after the 10,000 substitutions. 
Expressed in another way, Table 2 shows the number of dots in the white area in Fig.2. 
In Fig. 7, each circle in the complex plane shows the value of k to meet the requirement 
A(k) — B(k) — 1, and the total number of circle in each figure corresponds to the number 
in Table 2, obviously. From these figures and tables, it is understandable that the sign of 
m u is not important. Perhaps the reason is that m u is very small, and it is almost negligible 
in comparison with other masses. 



III. CONCLUSION AND DISCUSSION 

In conclusion, we have discussed the probability that the following model will be realized 
without fine tuning. The random numbers which become normal distributions have been 
substituted for each physical value at ji = m Z - And we have taken the RGE effect between 
\i = m z and Ax into consideration. In this way, the search for k which sets A(k) and 
B(k) to 1 simultaneously has been repeated 10,000 times. By this way, we have arrived at 
three conclusions: (1) The probability that the model will be realized without fine tuning is 
about 5% if we select the appropriate signs (14) or (15) of the masses. (2) This probability 
will increase if the signs of rrid, and m s are same. This gives the suggestion to the texture 
model. For example, a model with a texture (M d ) n = on the nearly diagonal basis of the 
up-quark Yukawa coupling M u is denied because these model leads to rrid/m s < 0. (3) From 
Fig.3-Fig.6, this probability will increase if we make m s somewhat larger or smaller than 
the present experiment value properly. 

In the present paper, we have demonstrate that the quark and charged lepton Yukawa 
coupling can be unified into only two matrices. However, we have not referred to the neutrino 
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masses and lepton flavor mixings. The neutrino Yukawa coupling is given by 

M° D = ^c u0 M° + C -^c ul M°, (3.1) 

CdO c dl 

Ml = c- R lM° D M^M° D T . (3.2) 



Concerning this problem, we have not been able to find the positive solutions within 3a 
which is written by the paper [5] for the present. However, since there are many possibilities 
for the neutrino mass generation mechanism, we are optimistic about this problem. 
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TABLES 





(m u , m c , m t ) 


(m d , m s , m b ) 




(m u , m c , m t ) 


(m d , m s , m b ) 


(0) 


(+ + +) 


(+ + +) 


(8) 


(+ + +) 


(+"+) 


(1) 


(- + +) 


(+ + +) 


(9) 


(- + +) 


(+"+) 


(2) 


(+-+) 


(+ + +) 


(10) 


(+-+) 


(+"+) 


(3) 


(--+) 


(+ + +) 


(11) 


(--+) 


(+-+) 


(4) 


(+ + +) 


(" + +) 


(12) 


(+ + +) 


(--+) 


(5) 


(- + +) 


(" + +) 


(13) 


(- + +) 


(--+) 


(6) 


(+-+) 


(" + +) 


(14) 


(+-+) 


(--+) 


(7) 


(--+) 


(" + +) 


(15) 


("-+) 


("-+) 



TABLE I. The combinations of the signs of (m M , m c , m t ) and (m^, m s , m^. The signs of the 
charged lepton are negligible in Eq.(2.4). 





sum 




sum 




sum 




sum 


(0) 


344 


(4) 


34 


(8) 


56 


(12) 


283 


(1) 


328 


(5) 


30 


(9) 


60 


(13) 


294 


(2) 


225 


(6) 


35 


(10) 


54 


(14) 


470 


(3) 


209 


(7) 


35 


(11) 


56 


(15) 


482 



TABLE II. The total number of the cases conforming to the requirements Min(i?(fi)) < 1 < 
Msx(B(k)) after the 10,000 substitutions. 



8 



FIGURES 
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FIG. 1. The relations in Eq.(2.4) on the complex plane of k. The solid line show A(k) = 1 and 
the dotted line B(k) = 1. This is an example which is given as follows: D® = diag(— 4.5116 * 1CT 6 , 
-0.0011789, 0.5028), D° d = diag(-5.6008 * 10~ 5 , -0.0007776, 0.038776), D° e = diag(1.8697 * 10~ 5 , 
0.0039461, 0.067375), 0? 2 = 0.22695, 9% 3 = 0.035057, 6^ = 0.0023936, and 5° = 1.3173. 
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FIG. 2. The maximum and minimum of B(k) on the line of A(k) = 1 in the case of (15) in 
Table 1. There are 10000 dots in all area, and 482 dots in the white area as tabulated in Table 2. 
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FIG. 3. The histograms show the distribution of data values conforming to the requirements 
Mm(B(n)) < 1 < Mslk(B(k)). Each number in parentheses show the signs of the mass in Table 
1. We bins the data values into 20 equally spaced containers, and show the number of elements in 
each container as a bar graph. The vertical solid and dotted lines show the center value and range 
of error in Eqs.(2.6) - (2.12), respectively. 
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FIG. 4. The histograms show the distribution of data values as Fig. 3. 
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FIG. 5. The histograms show the distribution of data values as Fig. 3. 
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FIG. 6. The histograms show the distribution of data values as Fig. 3. 
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FIG. 7. The distribution of k in the complex plane. Each circle shows the value of k to meet 
the requirement A(k) = B(k) = 1 
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